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We present a preliminary study of the neutral kaon mixing bag parameter Bk using two flavors of 
dynamical Wilson fermions. We determine the matrix element of the relevant AS = 2 operator by 
using both the conventional approach and the so called "non-subtraction method", and find that 
the latter leads to results with smaller uncertainties. After having implemented non-perturbative 
renormalization, we study the dependence of Bk on the see quark mass. At our relatively heavy 
values of quark masses (Mp/My — 0.60-^0.75) such a dependence is found to be negligible and 
the results, within the statistical accuracy, are consistent with a quenched determination. As a 
preliminary result for the renormalization group invariant parameter we quote Bk = 1 .02(25). 
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1. Preamble 

The parameter Bk(h), defined through 

(Fl^V) I K°) = ^F£m 2 K B K (LL) , (1.1) 

with in „ 

G^ 2 =%(1 -y 5 )^(l -75M (1.2) 

includes the long distance QCD effects in the evaluation of the indirect CP-violating parameter 
%. The theoretical uncertainty in the determination of B K represents at present one of the main 
sources of error in the unitarity triangle analysis [JlJ]. Within the quenched approximation the lattice 
estimates of Bk have reached an accuracy better than 10% [Q]-[|6|], so that the quenching effect, 
which has been estimated to be as high as ±15% [^], remains the primary source of systematic 
error to be investigated. Preliminary unquenched studies are presented in [Q] . 

In this talk we report on a exploratory calculation using Nf = 2 degenerate flavors of dynamical 
fermions. With respect to previous calculations, we simulate lighter sea quark masses on a finer 
lattice, but we work with unimproved Wilson fermions. We use the Wilson plaquette gauge action 
and the Wilson quark action at /3 = 5.8, which corresponds to a lattice spacing cC x ~ 3.2 GeV. The 
lattice volume is 24 3 x 48. The numerical simulation is performed by using the Hybrid Monte Carlo 



(HMC) algorithm [10]. A sample of 50 configurations have been generated at four different values 



of sea quark masses, for which the ratio of the pseudoscalar over vector meson masses lies in the 



range Mp/My z± 0.60^-0.75. Further details on the numerical simulation can be found in [1 1]. Our 
kaons are made with degenerate valence quarks and we considered four values of valence quark 
masses for each dynamical quark. 

To extract Bk, the following three steps are needed: 

• renormalization of the relevant AS = 2 operators; 

• extraction of the matrix elements from the appropriate correlation functions; 

• determination of Bk and extrapolation to the physical quark masses. Since we use Wilson 
fermions, which break chiral symmetry at 0{a), cautions are required in this latter step [g]. 
In particular, in addition to the finite mixing with dim = 6 operators, we have to subtract the 
effects of dim>l operators at finite lattice spacing. 

2. Renormalization and Matrix elements 

The first ingredient in the evaluation of Bk is the renormalization of the g A5 = 2 operator of 
eq. (fT^). From eq. ( |1.1[ ), we see that only the parity even component of Q AS=2 , namely Q\ = 
sy^d sj^d + sYfiYsd sy^y<,d, gives a non- vanishing contribution to Bk- 

In regularizations with exact chiral symmetry, such as continuum dimensional regularization 
or Ginsparg-Wilson fermions on the lattice, the operator Q\ renormalizes multiplicatively. For 
Wilson fermions, instead, the renormalization pattern is more involved, and can be expressed as 



Gi(/i) =Z vv+AA (ayi) 



5 



0i(a) + £A,(a)Gi(a) 

i=2 



(2.1) 
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Figure 1: The renormalization constant Zyy' +AA obtained with the RI-MOM method. A residual G((an) 2 ) 
dependence is observed. Our value of Zyy' +AA in the text is given by the intercept of a linear fit in (ajj.) 2 
(black line). 

Here Zw+aa(^I^) is the multiplicative renormalization constant, present also in formulations where 
chiral symmetry is preserved, while A2-s(a) are mixing coefficients peculiar for the Wilson regu- 
larization. The corresponding four-fermion operators are 



sd sd + sy$d sy$d 



so^ v d so^ v d. 



(2.2) 



g 2 (ju) = sy^d s-fd-sy^d sy^y 5 d , Q 3 (m) 
2 4 (a0 = sdsd-sy 5 dsy 5 d , Q 5 (n) 

In this study the renormalization constants Zw+aa^I^) and the mixing coefficients A2-s(a) have 



been computed non-perturbatively with the RI-MOM method [12]. The procedure is illustrated in 



details in [ ]13| ] for the quenched case. In the unquenched case, the situation is similar. The only 
additional constraint is that, in order to get these constants in a mass-independent renormalization 
scheme, a chiral extrapolation in both the valence and sea quark masses has to be performed. At 
this stage, the Goldstone pole effects have been non-perturbatively subtracted [O], even though 
their contributions turn out to be negligible. Moreover, in order to get rid of potential &{[a\i) 2 ) 
lattice artifacts, we have performed a linear fit in (aju) in the range [0.8, 2.5]. As an example, 
we show this procedure in fig. [T] for the renormalization group invariant (RGI) constant Zyy ! +AA . 
The scheme and scale dependent Zw+aa{^) and/or fix-(ju) are related to their RGI expressions 
through 



B K = [a s Gu)] 



2ft 



4tt 



Bk(h), 



(2.3) 



11 —2Nf/3. We use the RI-MOM scheme for which the NLO coefficient 
= 2. Our results at j3 =5.8 for the renormalization constant and mixing 
= 0.75(6), A 2 = -0.08(10), A 3 = -0.05(3), A 4 = 0.02(4), A 5 = 0.03(2) . 



where y = 4 and j3 = 
J = 2.83551 with N f = 
coefficients are Zyy ! +AA 

As known from quenched studies [12], the practical difficulty in the calculation of Bk with Wilson 
fermions is not only the non-perturbative evaluation of the subtraction constants A2_5(a) but also 
the necessity to perform it with a high level of accuracy, since the lattice bare matrix elements 
{Q2-5) are orders of magnitude larger than (<2i). Therefore, even though the subtraction constants 
are numerically small, the net effect of the subtraction is large. 

An alternative approach that allows one to compute the matrix element ( |1.1| ) without the neces- 



sity to perform the subtraction has been developed in [14]. In this procedure "without subtractions" 
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Figure 2: The ratios R\ and R2 computed with mesons at rest. 

one uses the hadronic chiral axial Ward identity to relate the matrix element of the operator Q\ 
to its parity violating counterpart, J2[ = sy^y 5 d sy^d + sy^d sy^y 5 d . The latter only renormalizes 
multiplicatively with the constant ZypXa\l), which we have calculated with the RI-MOM method 
determining, Zy^ 1 = 0.80(2). In this way, the problem of mixing with the other dimension-six 
operators is circumvented. The price to pay is that one has to compute a four-point correlation 
function where one pion operator is integrated over all lattice space-time coordinates. 
In practice, we consider the following two ratios of correlation functions 

I xy (P(x,0Gi(0)P(y,f / ))^( x P + ^ 



P pq 



(0 



*r q (o 



E x (P(x,0i J (0))^I y (P(y,f / )P(0))^ 
I xyz (P(x,0(m AW// n(z))J 1 (0)P(y, ?/ ))^P+^) 



(2.4) 



(2.5) 



I x (P(x,0P(0))^PI y (P(y,f/)P(0))^ 

for the method with and without subtraction respectively. Thanks to the Ward identity introduced 
in 0], Rf q and only differ by @(a) effects. At large times t f , T - 1 » 0, the ratios Rj(t) are 
both proportional to the desired matrix element (P°|<2 A5=2 (aO | P°). In our simulation, we have 
chosen tf = 14 and t in the range [23,37], whereas the momentum configurations {p,q} (in units 
of27T/La) are given by {(0,0,0), (0,0,0)}, {(0,0,0), (1,0,0)} and {(1,0,0), (0,0,0)}. An average 
over momentum configurations equivalent under hypercubic rotations is also performed. In fig. ||, 
we compare the results for Pi and R2 in the cases of the heaviest and the lightest quark masses 
{nival = nisea)- We notice that R2 suffers from rather smaller statistical fluctuations. 

3. Chiral behavior at finite lattice spacing 

At finite lattice spacing, the matrix elements (P°|2 AS=2 (jU) | P°) extracted from the ratios P,- 
are no longer proportional to Bk- The best approach to get rid of the &{a) contributions coming 
from the mixing with dirndl operators would be to perform a continuum extrapolation. In the 
present simulation, however, we have data only at a single value of the lattice spacing. Therefore, 
we rely on the approach proposed in [|fl. By calculating matrix elements of external kaons with 
non-zero momentum, we introduce an additional degree of freedom which allows us to partially 
remove the leading lattice artifacts. By writing 



(P°|Gi 



P°,q) 



ttm sea + Pm se m 2 P + -Fjm 2 p y msea , 



(3.1) 
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Figure 3: Values of the coefficients a and /3 as obtained from the fit to eq. (3.1 ) as function of mp(m sea ) 



one finds that the physical contribution to Bk, at each m sea , is given by the coefficient 7, while 
a and j3 parameterizes pure lattice artifacts. The results of the fit for the coefficients a and j3 are 
shown in fig. <j3j> . These coefficients turn to be sizable and their contribution in the extraction of the 
parameter 7 cannot be neglected. It is also interesting to notice that the effect of these coefficients 
is less relevant in the method without subtraction (blue points in fig. (||)). 

We finally plot in fig. |]5jf(«w,) = y sea , as a function of the square of the pseudoscalar me- 
son mass mp(m sea ). We observe that the method without-subtractions (red points) looks more 
promising, by suffering from smaller uncertainties. We also compare in the plot our results with 
those obtained by the UKQCD Collaboration [j|] by using 0(a) -improved Clover fermions but at 
a larger value of the lattice spacing. The two sets of results are very well compatible. In order to 
obtain the physical value of Bk , an extrapolation to mp(m sea ) = should be performed. Within the 
errors, however, we do not see any significant dependence on the sea quark mass and our results, at 
the simulated values of quark masses, are compatible with the quenched estimate (yellow band in 
fig. Bk = 0.97(9) obtained by performing a quenched simulation on a lattice with similar size 
and resolution. As a om a constant fit to 
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Figure 4: Values of Bk, as obtained from the two methods with and without subtraction, as a function of 
the pseudoscalar meson mass squared. The results obtained by the UKQCD Collaboration H and the band 
corresponding to the quenched estimate are also shown for comparison. 
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the points obtained with the method without subtraction, namely 

5^ = 1.02(25) (3.2) 
In order to reduce the uncertainties, we plan to include in the analysis more external momenta by 



exploiting the 6 -boundary condition method of ref. [JT5p . 
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